Ultrafast pump-probe dynamics in ZnSe-based semiconductor quantum-wells 
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Pump-probe experiments are used as a controllable way to investigate the properties of photoex- 
cited semiconductors, in particular, the absorption saturation. We present an experiment-theory 
comparison for ZnSe quantum wells, investigating the energy renormalization and bleaching of the 
excitonic resonances. Experiments were performed with spin-selective excitation and above-bandgap 
pumping. The model, based on the semiconductor Bloch equations in the screened Hartree-Fock 
approximation, takes various scattering processes into account phenomenologically. Comparing nu- 
merical results with available experimental data, we explain the experimental results and find that 
the electron spin-flip occurs on a time scale of 30 ps. 

PACS numbers: 190.0190, 190.7110, 190.5970, 320.0320, 320.7130. 
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I. INTRODUCTION 



From an experimental point of view, one can investi- 
gate the optical properties of semiconductors by exciting 
carriers (by means of optical pumping or carrier injec- 
tion) and measuring the absorption of a subsequent probe 
pulse. By comparison of this spectrum with the linear 
absorption spectrum, one obtains information on the in- 
fluence of the excitation on the absorption phenomenon 
and insight into the electronic and optical properties of 
the electron-hole plasma. Interpretation of experimental 
results is, however, nontrivial, given the substantial in- 
fluence of Coulomb and many-body effects, which give 
rise to a rich variety of broadening and energy renor- 
malizations. Moreover, the time evolution of the initial 
electronholc plasma makes the whole problem challeng- 
ing, both theoretically and numerically. In this paper we 
present a model that can describe the time evolution of 
the nonequilibrium electron-hole system but that is also 
simple enough to account for many dynamical processes 
that occur, using different polarizations of the pump and 
probe beams. 

Much previous theoretical work has focused on the 
study the absorption phenomenon in semiconductor 

auantum-wells in a quasi-equilibrium situation [l|, |2, l3|, 
, Q- Here, based on the available experimental data, 
we move beyond such a quasi-equilibrium situation. We 
include six dynamical processes that lead eventually to 
a thermal quasi-equilibrium in the electron-hole plasma: 
relaxation of the hot carriers' distributions toward Fermi- 
Dirac distributions, thermalization among the carrier 
gases, plasma cooling, carrier spin-flip, scattering be- 
tween the light- and heavy- hole bands, and recombi- 
nation (both radiative and nonradiative). A true mi- 
croscopic treatment accounting for all these many-body 
effects would be computationally prohibitive. Instead, 
we use a phenomenological approach to describe the 
time evolution of the hot electron-hole plasma. The ab- 
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sorption spectra are evaluated from the time-dependent 
semiconductor Bloch equations (SBE) in the screened 
Hartree-Fock approximation. In this work, as we focus 
on time scales over many picoseconds, wc do not need 
to consider the full coherent dynamics involving the non- 
linear scattering of pump light into the probe direction 
as a result of the nonlinear polarization interaction. The 
focus we choose allows and justifies the phenomenologi- 
cal treatment of the scattering processes for our qualita- 
tive analysis. Varying the delay between the pump and 
probe beams will allow us to obtain the time evolution of 
both the bleaching of the exciton peaks and the energy 
renormalizations, which, compared to experimental data 
will give an estimation of the time scale of the scattering 
processes mentioned above. Simultaneously, as we show 
below, the dynamics of the electron-hole plasma (density, 
temperature, plasma screening and distribution of each 
type of carrier in each spin state) can be monitored. 

The aim of this paper is to present our model for 
the time evolution of the electron-hole plasma created 
by spin-selective excitation in the absorption continuum 
and study its influence on absorption spectra. The ex- 
perimental setup and results are described in Section 2. 
In Section 3, we present our theoretical model for the 
time evolution of the electron-hole plasma, including the 
semiconductor Bloch equations that have to be solved 
numerically together with the rate equations used. We 
discuss our numerical results, comparing them with ex- 
perimental data, in Section 4. 



II. ULTRAFAST PUMP AND PROBE 
EXPERIMENTS 

A femtosecond laser system consisting of a Beam- 
lok argon- ion (Ar ) laser, a Tsunami mode- locked 
Ti:Sapphire laser, a Merlin Q-switched Nd:YLF laser, a 
Spitfire pulsed Ti:Sapphire regenerative amplifier and an 
ultrafast kilohertz optical parametric amplifier (OPA), 
was used for the generation of the ultrafast pump pulses. 
The Ar^ laser and the Merlin laser were the excitation 
sources for the Tsunami laser and the Spitfire amplifier 
respectively. The Tsunami output was fed to the Spitfire 



where it was temporally stretched, amplified, and finally 
temporally compressed. The Spitfire output provided the 
pump beam for the frequency conversion processes in the 
OPA. The overall system was capable of delivering a IkHz 
train of ~ 150 fs pulses and the wavelength was tuned 
at 459 nm (~ 2.69 cV). For the generation of the white 
light continuum, ~ 5% of the Spitfire output (A = 800 
nm) was focused on a 10-mm-thick quartz cuvette con- 
taining deionized water. The pump pulse was used to 
excite the semiconductor sample, which was mounted on 
the cryostat and cooled at 4K. The pump power, and 
therefore the carrier density, could be controlled by the 
use of a neutral density filter. The pump pulse power in- 
cident in the cryostat was 0.06 mW. The changes induced 
in the transmitted probe pulse energy were measured by 
an optical spectrum multichannel analyzer as a function 
of the time delay between the pump and the probe pulses. 
Using a glass microscope slide to monitor its stability we 
selected a small portion of the white-light continuum be- 
fore it fell onto the sample. The spot size radius of the 
probe beam was 190 /zm, and it was considerably smaller 
than that of the pump in order to probe a region of uni- 
form photoexcited density. Both pump and probe beams 
were circularly polarized and independently controllable 
by A/4 plates. Opposite circular (OCP) and same circu- 
lar polarization (SCP) configurations were employed. We 
provided the time resolution by delaying the white light 
continuum pulses relative to the pump pulses. The exper- 
imental work was performed on a ZnSe/ZnCdSe multiple 
quantum well structure of twenty 4nm-wide wells grown 
by molecular beam epitaxy on GaAs substrate. The 20% 
Cd content in the wells produces light-hole-heavy-hole 
exciton splitting of more than 30 meV. 

The absorption spectra. Fig. 1, show that at early 
times, both the heavy-hole and light-hole exciton peaks 
are bleached but not shifted much. The broadening that 
is due to the interaction-induced dephasing plays an im- 
portant role in both SCP and OCP, but for SCP the 
heavy-hole exciton peak is more bleached because of the 
Pauli blocking effect that reduces the oscillator strength. 
In contrast, the light hole exciton peak is more bleached 
in OCP. A detailed description and an interpretation of 
these experimental results are given below with our nu- 
merical analysis. 

Increasing the delay between the pump and the probe 
beams suggests that many dynamical processes occur in 
the electron-hole plasma and change the shape of the ab- 
sorption spectra. The dynamics of the absorption spec- 
tra is shown in Fig. 2, where both the bleaching (Aa/a) 
and the energy shift of the heavy-hole exciton peak are 
given as functions of the delay between the pump and 
the probe. Experimental data show an overall decay of 
the exciton peak bleaching as well as a convergence of 
the OCP and the SCP curves. They also show an initial 
blueshift at early times and an energy shift that brings 
the resonances back to the linear spectrum exciton reso- 
nance. The energy shift exhibits the same type of behav- 
ior as the exciton bleaching: The OCP and SCP curves 
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FIG. 1: Measured absorption spectra. Comparison of OCP 
SCP absorption. 
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FIG. 2: Measured absorption spectra dynamics. Comparison 
of OCP and SCP heavy-hole exciton peaks' bleaching and 
shift. 



III. THEORETICAL MODEL 

To describe and explain what we observe, we con- 
structed a theoretical model that describes the time evo- 
lution of the electron- hole plasma and its influence on the 
absorption spectra. 



A. Polarization dynamics 

Inasmuch as the heavy-hole-light-hole band splitting 
is AE'cs = 30 meV, we neglect the heavy-hole (hh) 
and light-hole (Ih) coupling. The interband polarization 
equations ^ are: 



- *KfcW + <fe(i) - 1) fi^R,fcW ,(1) 
where the Rabi energies, huj^ j., are given by: 

for A = hh, Ih. The energies Bi^k are the renormahzed 
energies evaluated in the static plasmon-pole approxima- 
tion including the contribution of the pair continuum [l| : 



where e is the dielectric constant and P{lo) and £{iu) 
are the Fourier transform of the polarization function 
P{t) = dcv^kPk{t) and of the electric field £{t). It is 
a complex function whose imaginary part is proportional 
to the absorption: a(a;) oc Imx(a;) [1|]. 

To solve Eq. ([T]) one needs knowledge of the carrier 
distributions to evaluate the phase space filling factor 
and the plasma temperature and density to calculate the 
screened Coulomb potential energy V^ entering the def- 
inition of the Rabi frequency, Eq. ^. We neglect co- 
herent polarization nonlinearities, as we are considering 
above bandgap pumping and time scales longer than the 
dephasing time. 



ei,k — ii,k + '^cxcA{k) + - Ai?cH 



e, hh and Ih, (3) 



with Sexc,i the screened exchange self-energy, AE'ch the 
Coulomb hole energy Q, and n^^(t) are the occupancy 
of the carrier of type c with spin a at the time t. 

The temporal envelope of the probe field, £{t), is as- 
sumed to be Gaussian and the optical suceptibility x{^) 
is defined as: 



X(w) = 



e£{uj)'' 



(4) 



B. Evolution of the carriers distributions 



In our model the time evolution of the distribution 
function, n"^ ^i^), of a carrier c with a spin a takes sev- 
eral dynamical processes into account. Relaxation of the 
carrier distributions, carrier spin-flip, recombination and 
light hole scattering yield the following system of six cou- 
pled differential equations {cr = t or cr = ],, indicating 
the two spin-states): 
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(5) 
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These equations are numerically solved to yield the 
time dependence of the carrier distributions ^cfe(^) 
whose values are 0-1. The various terms that enter into 
the above system are now described. 

Intraband scattering (rcq terms): The optical 
pumping creates a population of hot carriers. One of 
the fastest processes (i.e. subpicosecond time scale Q) 
that occurs in each band is the rapid equilibration of 
these carriers: due to carrier-carrier scattering the initial 
hot carrier distributions evolve towards Fermi-Dirac 
quasi-equilibrium distributions. Extensive work on 



this specific topic involving the quantum Boltzmann 
equation can be found in the literature [^, llfl llll Il3 |- 
Here we use a phenomenological approach to describe 
the time evolution of the hot carrier distributions which 
is characterized by a relaxation time Toq 
with intraband scattering. Thus, the quantities n^ ^. 
in Eq. ([5]) are Fermi-Dirac distributions describing the 
quasi-equilibrium for each spin-polarized subsystem, 
which has the same carrier density and energy as the 
nonequilibrium distribution. Note that the intraband 
scattering does not change the carrier densities nor the 



associated 

eq,(T 



total kinetic energy. 



plasma density. 



Carrier thermalization (rthcim terms): A process 
that also influences the time evolution of the carrier 
distributions is the thermalization process among carri- 
ers of different types. The scattering between electrons 
and heavy- and light-holes is a process that drives the 
initial carrier temperatures, T^, to a common quasi- 
equilibrium temperature, Toq, which can be different 
from the lattice temperature, Tiat- To evaluate Toq, 
one needs to calculate the total plasma energy Etot 
and then compute the corresponding temperature Toq 
assuming a Fermi-Dirac distribution. To account for 
the thermalization process, we suppose that the time 
evolution of the distribution functions is characterized 
by a phenomenological time rthcrm- Thus, the quantities 
in Eq. (O are Fermi-Dirac distributions that 



therm. (T 



describe the quasi-equilibrium for each spin-polarized 
subsystem; they have the same carrier density but an 
energy corresponding to the common quasi-equilibrium 
temperature Teq. 

Carrier spin-flip {tsi terms): The spin- flip is a 
process by which the spin orientation of a carrier is 
reversed. Models that describe such a phenonemon 
have been proposed^ e.g. Elliot- Yaffet [ij, [lj| and 



D'Yakonov-Perel [15|, but the detailed mechanisms 
responsible for the spin-flip process are not yet well 
understood despite extensive studies [1^, [l3, [3 ■ In this 
work, we only consider the spin-flip phenomenologically 
with the associated characteristic time Tsi that is of 
interest to us. The spin- flip introduces a coupling 
between the spin-states a and a' for a given type of 
carrier c, which considerably complicates the solution of 
Eq. ([5]) . The value of Tgf is poorly known and one of our 
aims using this model will be to extract this value from 
the experimental data. 

Recombination (xiad and r^r terms): We distin- 
guish here between radiative and nonradiative recombi- 
nations. Radiative recombination is a process that oc- 
curs on relatively long time scale compared to the others 
described above (T^ad = 1.6 ns in ZnSe). The total ob- 
served luminescence intensity /pi(t) is directly linked to 
the distribution functions of these carriers: 



^p'iW«E 



Mcvr <Jt) nlJt) , 



(6) 



for the heavy and light hole optical transitions, where 
A = Ih or hh, and dcv is the dipole matrix element. 

For a detailed study of radiative recombination and 
spontaneous emssion rate, see Ref. [la|. Detailed 
calculations with which to evaluate Trad can be found in 
Ref. |20!| . The term Tnr in the right hand side of Eq. ^ 
accounts for the nonradiative recombination which is a 
faster process than the radiative recombination at high 



Heavy hole and light hole scattering (rjh terms): 
Away from zone center, both heavy and light holes are 
mixtures of the bulk valence band states. This mix- 
ture enhances the intersubband scattering between the 
heavy- and light- hole bands. Many processes that involve 
other quasi-particlcs (heavy and light holes, phonons, ex- 
citons . . . ) can facilitate this type of scattering, and it 
is nontrivial to assess the relative importance of each of 
the processes. To avoid such complication we make a 
simple approximation, assuming that this scattering is 
spin independent, e.g., that |3/2, l/2)ih scatters equally 
into |3/2,3/2)hh- When the quasi-equilibrium between 
the heavy-hole and the light-hole bands is reached, their 
chemical potentials have to satisfy the relation /zih — 
^hh — AE'cs because of the band splitting that is due 
to confinement and strain. Solving this equation at low 
temperature and for typical electronhole plasma densi- 
ties, i.e. of the order of lO^^cm^^ or more, yields a neg- 
ligible light-hole density: N\\i/N « 0. Hence we include 
a simple decay characterized by the time Tih in Eqs. (5) 
to model the light-hole scattering. 



C. Plasma cooling 

Intraband scattering dominates the fast carrier distri- 
bution relaxation but is not a process that dissipates en- 
ergy. Thus the initial electronhole plasma temperature is 
determined by the kinetic energy of the nonequilibrium 
distribution created by the femtosecond pump pulse. De- 
pending on the energy of the excitation pulse, the effec- 
tive plasma temperature can be well above the lattice 
temperature. The most important source of energy dis- 
sipation is the coupling of the electronic system with the 
lattice. The plasma cooling can be treated by solution 
of the quantum Boltzmann equation that describes the 
carrier-phonon scattering [ll|, IT4I, but we restrict our- 
selves to a phenomenological approach. Hence, the loss 
of carriers' kinetic energy obeys a rate equation: 



d-£^tot,e 

dt 



pCq,(T 



^tot.c 



(7) 



where i?tot"c is the total energy per unit area calculated 
from the quasi-equilibrium distribution i\^^ at the lat- 
tice temperature. Eq. ([7]) has to be solved numerically. 

To evaluate the effective carrier temperatures at each 
point in time, we use the time evolution of the nonequi- 
librium energies E^^^,.{t). We calculate first the quasi- 
equilibrium energies as an explicit function of tem- 
perature E^'i'" (T;^); then, equating El^^^" {T^) and the 
nonequilibrium energies E^^^ ^(t) gives the effective tem- 
perature T^ at the time t: 
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1 



,(8) 



where /3^ = l/fcsTf . 

Knowledge of the distributions n'^^i^) and tempera- 
tures T^ allows us to evaluate the time evolution of the 
Pauli blocking factor that enters into the equations of 
motion for the interband polarizations, Eq.([T]). 



D. Time evolution of the plasma screening 



E. Initial conditions for the electron/hole plasma 

We now turn to the three pump-probe polarization 
configurations and the initial carrier distributions they 
produce as the starting point for the solution of Eqs. ([5]): 

1. OCPs, for which the pump and the probe beams are 
both circularly polarized but in opposite senses, 

2. SCPs, for which the pump and the probe beams 
are cicularly polarized in the same sens, 

3. SLPs, for which the pump and the probe beams are 
both linearly senses. 



As the plasma temperature and the carrier densities 
evolve, the plasma screening changes. In this paper we 
are concerned with the time evolution of the electron-hole 
plasma and we need to find a simple way to evaluate it 
in a nonequilibrium situation. This Lindhard formula, 
which is valid for both equilibrium and nonequilibrium 
situations [l], is our starting point: 



(l.) 



y.ET. 



h{ui — iS) + e 



(9) 



c^k—q 



^c,k 



For ease is notation we do not explicitly denote the 
time dependence of the various physical quantities de- 
fined above. We have to simplify Eq. ^ as it is not 
practical to use that equation for numerical purposes be- 
cause of its continuum of poles. We choose to work in the 
long wavelength limit q — *■ [1|. Using a nonequilibrium 
distribution function we find: 



e(g->0,w) = l 



N° 



^1 ,„2 Z^ ^„ 



(10) 



which shows that the time dependence of the distribu- 



tion functions 



'■cM 



due to the relaxation process does 



not affect the expression of e((7 — > 0, w): it has the same 
form as in quasi-equilibrium calculations. So, to treat the 
nonequilibrium problem for the screening, we can calcu- 
late the screening from the quasi-equilibrium formulas 
that one can find in Refs. [l|, l2l|. The time dependence 
of t{q -^ 0,oj) is contained in N^ . To obtain analyti- 
cal results for the plasma screening we make use of the 
static plasmon-pole approximation [ij, and the screened 
Coulomb potential given by: 



1 



(11) 



Note that our definition of the bare 2D Coulomb po- 
tential in a quantum well includes the form factor fq to 
account for the finite well width w [3|. 
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FIG. 3: Selection rules for optical transitions achieved with 
polarized pump light 

Depending on the pump-probe polarization config- 
uration, we generate different couplings between the 
ground state and the various spin-states |J, rrij). In our 
analysis we include the two-fold degenerate conduction 
band |l/2, l/2)e and |l/2, -l/2)e, the two-fold degener- 
ate heavy hole band \3/2,3/2)iih and \3/2,—3/2)hh and 
the two-fold degenerate light hole band |3/2, l/2)ih and 
3/2, — l/2);/i. The selection rules for zinc blende semi- 
conductors are used [231 and the relative populations gen- 
erated by optical pumping in the continuum are as de- 
picted in Fig. 3. The ratio of 3 between the populations 
created from the heavy hole and light hole transitions 
comes from the ratio of the dipole matrix elements which 
describe the relative strengths of these optical transitions 

In principle, the carrier dynamics during the pump 
process could be included in our model. However, since 
Tcq ~ 100 fs and we are not interested in the early co- 
herent regime, we simply assume an initial carrier den- 
sity with nonthermal Gaussian distributions reflecting 
the spectral width and the location of the pump pulse. 
This assumption has no significant impact on our nu- 
merical results as Tcq is of the order of 100 fs and saves 
considerable computational effort. In the SLP situation, 
electrons with a given spin a are created from both the 
light- and heavy-hole transitions. So, the initial distri- 
bution for the spin-polarized electron gases in the SLP 
is the sum of two Gaussian distributions, each centered 
appropriately. 



IV. COMPARISON BETWEEN NUMERICAL 
AND EXPERIMENTAL RESULTS 

To present and discuss the dynamics of the absorp- 
tion spectra, it is useful first to describe the time evolu- 
tion of the electron-hole plasma, which is influenced by 
the dynamical processes mentioned above. We are con- 
cerned with the interplay among the various dynamical 
processes included in our model and their influence on 
the time evolution of the carrier density, plasma energy 
and temperature, situations, namely, that when excita- 
tions have been performed by circularly polarized light 
and that when they have been performed by linearly po- 
larized light. The case of circular polarization (CP) can 
describe both OCP and SCP situations: The dynamics 
of the carrier gas is the same; only the interaction with 
the probe field is different. 

The lattice temperature is taken to be Tiat = 77 K (we 
find in materials with strong Coulomb effects such as 
ZnSe, that strongly degenerate cases at 4 K are numer- 
ically prohibitive owing to the large number of /c-states 
required near the Fermi energy) and the initial plasma 
density A^ = 3 x lO^^cm"^. The effective masses we use 
are rUc = rTijh — 0.15?no and mhh = O.Gttiq, where toq 
is the mass of the free electron. The dielectric constant 
is e = 8.8 and the bandgap Eg — 2.66 eV. The phe- 
nomenological parameters entering Eqs. ([5]) are: r^q — 
0.1 ps, Tcooi = 1 ps, Tthcrm = 1 ps. Thcsc characteristic 
times only give an order of magnitude and are taken from 
Ref. ^]. The radiative recombination time Trad = 1-6 ns 
is calculated for ZnSe parameters. We assume the non- 
radiative recombination to be density independent with 
Tni- — 30 ps taken from experiments, by fitting numerical 
results to experimental results (see Figs. 2 and 8). Inas- 
much as we are pumping in the continuum, the excited 
valence band states are of a mixed heavy-hole-light-hole 
character and hence scatter efficiently. We chose rih = 0.5 
ps for the light hole density decay. The spin-flip times, 
Tgf = 30 ps, have been chosen also comparing numerical 
results with experimental data. Various values for Tn,. 
and Tsf have been tried numerically and we estimate the 
margin of error for the parameters to be ±5 ps. 



Time evolution of the carrier densities 
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FIG. 4: Time evolution of the spin-polarized gas popula- 
tions (lO^^cm" ) created by spin-selective excitation. Note 
the short time scale for the light hole density decay. 



That is the reason why the initially less populated 
electron gas exhibits a maximum at early times (10 
ps). Note that for SLP, the spin-flip process has no 
effect on the dynamics of the electron densities. Hence, 
the only contributor to the density decay in the SLP 
situation is recombination. Also note that the radiative 
recombination is too slow (1.6 ns) to have much influence 
on the fast population dynamics. 



Heavy holes: Similar comments apply for the heavy 
hole population dynamics as above with t^^ = 30 ps. 
However, the fast intersubband scattering that drives 
the light holes into the heavy hole band on a time scale 
given by rih = 0.5 ps has to be considered here. That 
explains the fast initial rise of both heavy hole densities, 
before they start decaying. Note that the initial spin 
down heavy hole density is zero, unlike the electrons 
whose populations in both spin states is non-zero from 
the beginning because of optical transitions from both 
the heavy- and light-hole bands. 



First we study the time evolution of the carrier den- 
sities obtained by solution of Eqs. ([5]) and shown in Fig. 4. 

Electrons: The processes responsible for the change 
in population are the electron spin-flip and the recombi- 
nation. The population of spin down electrons decreases 
while the population of electrons with opposite spin 
increases. The carrier spin-flip is a process that tends 
to create equal spin populations over a time scale given 
by r°f = 30 ps. Then, the recombination which occurs 
on the same time scale given by Tnr = 30 ps becomes 
dominant and drives both populations towards zero. 



Light holes: The time evolution of the light hole pop- 
ulations can also be described similarly as the electron 
populations with r^f = 30 ps. However, as in the case 
of the heavy holes, one has to consider the intersubband 
scattering. As it is a very fast process, the population of 
light holes decreases on a very short time scale given by 
Tih = 0.5 ps. After one picosecond the light hole density 
is negligible (but non-zero as heavy holes still scatter into 
the light hole band). Note that in the case of the spin 
down light hole gas, the population always remains very 
low: the spin-flip time is comparatively too long to make 
any significant change. 



B. Time evolution of the plasma energy 



C. Time evolution of the plasma temperature 



The time evolution of the plasma energy is displayed 
on Fig. 5. 



The effective temperatures of each carrier gas, shown 
in Fig. 6, are computed solving Eq. ([5]). 
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FIG. 5: Time evolution of the spin-polarized gases energy loss 
(eV//ini ). Note the very short time scale for the light hole 
gases' energy loss. 



FIG. 6: Time evolution of the spin-polarized gases' temper- 
atures (K). Note the very short time scale for the light hole 
temperature evolution. 



Electrons: The initial total energy of the spin down 
electron gas in the CP case corresponds to a temperature 
above the lattice temperature and hence is a decreasing 
function of time because of the cooling process that 
occurs on a time scale given by t^q^j = 1 ps. One can 
observe that the energy of the electron gas with opposite 
spin is an increasing function of time before cooling down 
to the lattice energy. It is due to the thermalization 
process between carriers of different types and also to 
the fact that the density increases over the same amount 
of time because of the spin-flip , and hence increases the 
total energy. In the LP case, we only observe a cooling 
as the spin-flip has no influence on the time evolution of 
the densities. 

Heavy holes: The initial rise of the carrier densities 
for both heavy-hole gases as a result of light-hole 
scattering influences the behavior of the energy, which 
is an increasing function of time at early times. The 
thermalization process is also responsible for this rise, as 
the initial hole gases' energies are below the spin-down 
electron gas energy. In the case of the spin-down 
heavy-hole gas the increase occurs on a much longer 
time scale than for the opposite spin heavy-hole gas. 
This is the result of the spinflip process that increases 
the population for the spin-up heavy-hole gas in CP. 

Light holes: The initial rise of the energy for the 
light hole gases is only due to the thermalization process. 
As the light hole density decay is fast, it only takes a 
few picoseconds for the light hole gas energy to become 
negligible. The population of the spin down light hole 
gas remains very low; so does the energy. 



Electrons: The time evolution of the electron 
temperatures follows exactly the time evolution of 
the electron gas energies, but it converges fast toward 
the lattice temperature, whereas the electron energies 
keep decreasing. This behavior is due to the fact that 
the electron densities also decrease by recombination, 
keeping the average electron kinetic energy constant. 

Heavy holes; The heavy-hole gas energies increase at 
early times; so do the densities, because of the light-hole 
scattering. In CP, for the spin-down heavy-hole gas, 
despite this increase in energy, the average heavy-hole 
kinetic energy decreases. Hence the effective tempera- 
ture decreases. Then, because of the thermalization with 
the electrons and the lattice, the temperature starts 
increasing. For the spin-up heavy-hole gas we observe a 
monotonic increase toward the lattice temperature as it 
gains energy from the thermalization processes with the 
other carriers and the lattice. However, as the densities 
starts decreasing, the energy of the heavy holes also 
starts decreasing. But, as in the case of electrons, there 
is compensation between the energy loss and the density 
decay that makes the average heavy-hole kinetic energy 
constant when it reaches the lattice temperature. 

Light holes: The light-hole density decay is so fast 
compared to the plasma cooling that the average light- 
hole kinetic energy keeps increasing and the light-hole 
gases temperatures go beyond the lattice temperature. 
In fact, the light-hole gases have no time to reach quasi- 
equilibrium with the lattice. We stop calculating the 
light-hole gases' effective temperatures when their popu- 
lation is small enough and their temperature high enough 



to have no influence on the absorption spectra (after 1.5 
ps). 



D. Time evolution of the absorption spectra 

In this section we present numerical solution of Eqs. ([1]) 
and ([5]) and discuss the behavior of the excitonic peaks' 
bleaching as well as their energy shift. First, we compare 
the calculated absorption spectra in Fig. 7 with exper- 
imental data on Fig. 1 for a given delay. The optical 
pumping is set 30 meV above the band edge, in the con- 
tinuum, thus creating an initial unbound electron-hole 
plasma. 
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FIG. 7: Calculated absorption spectra. Comparison of OCP, 
SLP and SCP absorption. 



The experimental and numerical spectra look similar, 
but one can observe a significant redshift of the OCP exci- 
ton peak that is not observed on experimental data. This 
artificial redshift is mainly due to the screened Hartree- 
Fock approximation [23] and also to the screening model 
leading to an overestimation of the bandgap renormaliza- 
tion Q. The heavy-hole exciton peak is more bleached 
in SCP than it is in OCP, and it is more blue shifted. 
This is due to the phase space filling cff'ect that is more 
important in the SCP configuration. For the light-hole 
exciton, as in the experiment, we observe the opposite 
phenomenon: the OCP light-hole exciton peak is more 
bleached and blue shifted than the SCP light hole exciton 
peak. This is because the spin-up electrons excited with 
the (7^- polarized pump, from the heavy hole transition, 
occupy states that would be created from a light-hole 
transition with the CT+-polarized pump. So, in the OCP 
situation, the reduction of the oscillator strength owing 
to filling of the phase space when one probes the light- 
hole transition is more important than it is for SCP. The 
SLP configuration can be seen as an intermediate case 
between OCP and SCP. 



E. Heavy hole exciton peak dynamics 

As we mentioned above, we are concerned with 
the time evolution of the absorption spectra. So we 
computed both the bleaching and the energy shift of 
the exciton resonances as functions of the time delay 
between the pump and the probe. The numerical results 
in Fig. 8 show a decay of the heavy-hole exciton peak 
bleachings and energy shifts for OCP, SLP and SCP. The 
three curves converge after a few tens of picoseconds, 
because of spin-flip and recombination. 
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FIG. 8: Calculated absorption spectra dynamics. Comparison 
of OCP, SLP and SCP heavy-hole exciton peaks' bleaching 
and shift. 

OCP: The initial state of the electronhole plasma 
given in Fig. 3 shows that the only contribution to heavy- 
hole exciton bleaching is due to plasma screening, which 
lowers the Coulomb enhancement. The rapid equilibra- 
tion of the carrier distributions, together with plasma 
cooling, contributes to the increase of bleaching at early 
times, i.e., less than 5 ps. Although in this case the 
carrier spin-flip occurs on the same time scale as the re- 
combination, the latter process is dominant, and, with 
decreasing plasma density plasma screening and phase- 
space filling become less important. Thus the amount 
of bleaching decreases. In terms of energy shift, the ini- 
tial heavy-hole exciton peak is redshifted. The bandgap 
renormalization is not strong enough to compensate for 
the exciton binding energy, which remains large because 
of the absence of the Pauli blocking effect. The redshift 
becomes even more important as the plasma screening is 
increased, because of the fast light-hole scattering. Then, 
because of the carrier spin flip and the recombination 
process, the heavy-hole exciton peak shifts toward the 
blue region and saturates on a longer time scale (from 30 
ps). This behavior is qualitatively different from what 
we observe in the experiments because of the initial large 
redshift: If the initial heavy- hole exciton peak were not 
redshifted it would enter the blue region because of the 
increasing phase-space filling effect owing to the spin flip 



and to plasma cooling; then, because of recombination, 
we would observe a shift toward the red region, which 
would explain the presence of a maximum. 

SCP: As well as plasma screening, phase-space fill- 
ing contributes to bleaching of the heavy-hole exciton 
peak. Because of thermalization and plasma cooling, we 
observe a slight increase for short (< 2ps) pump-probe 
delays, but the carrier spin flip decreases the amount of 
bleaching together with the recombination process. In 
SCP both processes contribute to decreasing the Pauli 
blocking filling effect. The initial SCP heavy- hole exci- 
ton peak is blueshifted, even though the bandgap renor- 
malization is more important than it is for OCP (the 
exchange term is nonzero) . This is so because the phase- 
space filling factor decreases the exciton binding energy 
significantly in the SCP case. Because of the fast re- 
laxation of the distributions and the plasma cooling, the 
heavy-hole exciton blueshift is even more important at 
early times (< 2ps); then the spin flip and recombina- 
tion make the heavy-hole exciton peak shift toward the 
red region. In the SCP situation, carrier spinflip and re- 
combination have the same effects on the exciton peak 
dynamics, whereas in the OCP case they tend to cancel 
each other out. 

SLP: As for SCP, plasma screening and phase-space 
filling contribute to heavy-hole exciton peak bleaching. 
There is an initial small increase that is due to thermal- 
ization and plasma cooling, but after a few picoseconds 
the recombination process starts effectively to influence 
the bleaching behavior. In SLP the spin-flip process plays 
no role in the dynamics of the bleaching. The behavior 
of the energy shift for SLP is similar to but less dramatic 
than the OCP behavior, as it is not influenced by the 
spin-flip process. 
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FIG. 9: Calculated absorption spectra dynamics. Comparison 
of OCP, SLP and SCP light-hole exciton peaks' bleaching and 
shift. 



exciton bleaching is characterized by a simple decay in 
OCP and SCP, whereas our calculations show the pres- 
ence of a maximum for short pumpprobe delay that is 
due to increased Pauli blocking because of plasma cool- 
ing. If the cooling characteristic time is increased, the 
numerical results show a plateaulike behavior for pump- 
probe delays of as much as several picoseconds before 
the recombination becomes dominant. The early (< 5 
ps) behavior shown in Fig. 8 for OCP is enhanced, and 
the SCP bleaching dynamics is qualitatively the same as 
that of OCP. 



DISCUSSION AND CONCLUSION 



F. Light hole exciton pealc dynamics 

In light hole transitions the contribution of the light 
hole population becomes negligible quickly as the light 
hole density decays on a vey short time scale, Tih = 0.5 
ps. The behavior of the light-hole exciton peak bleaching 
and shift is dominated by the electrons. As we mentioned 
above, the bleaching here is more important for OCP 
than for SCP. Because of the rapid light-hole density 
decay, the bleaching decreases at very early times (less 
than 1 ps), but it increases shortly after that because of 
electron thermalization and plasma cooling, as shown in 
Fig. 9. Then, after a few picoseconds, the recombination 
process dominates the behavior of the bleaching, which 
decreases monotonically with increasing pumpprobe de- 
lay. The three curves converge in ^^30 ps. As far as the 
energy shift is concerned, the blueshift quickly becomes 
small because of the fast light-hole density decay. 

Calculations made with longer characteristic times for 
the plasma thermalization and cooling (> 5 ps) led to re- 
sults that are qualitatively quite different from exprimen- 
tal data: In Fig. 2 the time evolution of the heavy-hole 



We have presented a model that describes the time evo- 
lution of an initial hot electron-hole plasma and its influ- 
ence on absorption spectra. We considered spin-selective 
excitation and included various dynamical process in our 
model. Except for the energy shift of the OCP heavy- 
hole exciton, the numerical results are in good qualita- 
tive agreement with the experiment. According to our 
calculations, the electron spin populations equilibrate on 
a much longer time scale (30 ps) than the thermalization 
of the electron-hole plasma (1 ps). 

This spin-flip time is consistent with the experimental 
values for the II— VI quantum wells discussed in Ref. [2J] . 
One would need more-extensive experimental data to ob- 
tain more conclusive results for the heavy-and light-hole 
spin-flip times. We made additional numerical calcula- 
tions to explore various spin-flip time regimes. The re- 
sults show that, when the electron spin- flip time is short 
(1 ps), the OCP and SCP curves converge faster than 
when it is long, even with a short heavy-hole spin-flip 
time. This means that the electrons, which are lighter 
than the heavy holes, always dominate phase-space fill- 
ing, even if the initial electron gas temperature is much 
higher than the heavy hole gas's temperature (Fig. 6). 
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This implies that our results are rather insensitive to the 
specific value of the heavy-hole spin flip time. The light 
holes do not play a siginifcant role in phase-space filling, 
as they decay on a very short time scale (0.5 ps). Plasma 
cooling enhances plasma screening, and together with the 
relaxation of the distributions it also enhances the Pauli 
blocking effect. We can observe the influence of these 
fast processes on the early behavior of the exciton peaks' 
bleaching and energy shift. The spin- flip process leads to 
opposite qualitative behavior of the bleaching dynamics, 
depending on the initial polarization configurations. As 
shown above, exciton peak bleaching and shift can either 
increase or decrease because of the spin-flip process, de- 
pending on the type of exciton (heavy or light) and on 
the polarization (OCP or SCP). The radiative recombi- 
nation occurs on a time scale that is too large ('^1.6 ns in 
ZnSe) to have any effect on the fast bleaching and energy 
shift dynamics (below 100 ps). However, the nonradia- 
tive recombination is fast enough that an overall decay 
can be observed on a time scale shorter than 100 ps. 

As far as the energy shift is concerned, the model that 
we used for this study seems not to be good enough to 
describe qualitatively the OCP and SCP energy shifts. 
We make two comments on this problem: First, we have 
restricted ourselves to the screened Hartree-Fock level. 



neglecting additional correlation terms to account for 
broadening, which is described here by phenomenological 
damping. Inclusion of such terms has been shown to cor- 
rect spurious shifts introduced by the screened Hartree- 
Fock approximation [23,] in SLP. However, the static 
plasmon-pole approximation that we and many others 
have used is constructed in just such a way as to produce 
a constant exciton energy at the screened Hartree-Fock 
level in SLP. We do this by adding a term to the effective 
plasmon frequency to account for the effects of the pair 
continuum. As can be seen from Fig. 7, there is effec- 
tively no shift of SLP in our calculations; we also checked 
our numerics against published GaAs results. Making 
the natural extension to OCP and SCP, however, leads 
to extraneous shifts, even when we include this pair con- 
tinuum term. Our conclusion is therefore that the static 
plasmon-pole approximation is inadequate for the treat- 
ment of spin-polarized populations. Of course one could 
alter the strength of the pair-continuum term in an ad 
hoc fashion to correct for this, but that correction would 
accomplish little. Within our model it is possible to con- 
trol the values of the phenomenological parameters to 
obtain further insight into the complex interplay among 
various dynamical processes as more experimental data 
become available. 
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